Commuting quantum transfer matrix approach to intrinsic Fermion system: 
Correlation length of a spinless Fermion model 
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The quantum transfer matrix (QTM) approach to integrable lattice Fermion systems is presented. 
As a simple case we treat the spinless Fermion model with repulsive interaction in critical regime. 
We derive a set of non-linear integral equations which characterize the free energy and the correlation 
length of {c'-Ci) for arbitrary particle density at any finite temperatures. The correlation length is 
determined by solving the integral equations numerically. Especially in low temperature limit this 
result agrees with the prediction from conformal field theory (CFT) with high accuracy. 
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I. INTRODUCTION 

Exact evaluations of physical quantities at finite tem- 
peratures pose serious difficulties even for integrable 
models. One has to go much beyond mere diagonaliza- 
tion of a Hamiltonian; summation over the eigenspectra 
must be performed. p, p, 

The string hypothesisETQ brought the first break- 
through and success. It yields a systematic way to eval- 
uate several bulk quantities including specific heats, sus- 
ceptibilities and so on. 

More recently, the quantum transfer matrix (QTM) 
method has been proposed to overcome some difficult 
to which the standard approach is not applicable.' 
One reduces the original problem to finding the largest 
eigenvalue of the QTM which acts on a fictitious sys- 
tem of size N (referred to as the Trotter number), which 
should be sent N ^ oo (Trotter limit). As this procedure 
is sometimes difficult, we integrate its procedure with an- 
other ingredient, the integrable structure of the underly- 
ing model. It allows for introduction of the commuting 
QTMs which are labeled by complex parameter xS3 A 
set of auxiliary functions, including the QTM itself, sat- 
isfy certain functional relations. We shall choose these 
functions such that they have a nice analytical property 
called ANZC (Analytic, NonZero, and Constant asymp- 
totics, see Sec. Ill) in a certain strip on the complex x 
plane. This admits the transformation of the functional 
relations into a closed set of the integral equations. For 
all cases known up to now, the Trotter limit N -^ oo can 
be taken analytically in the integral equations. We thus 
have seen a remarkable reduction from the problem of 
combinatorics (summation over the eigenspectra) to the 
study of analytic structures of suitably chosen auxiliary 
functions. 

This novel scenario ]jia£|-been applied to many mod- 
els of physical interest JiiTcZl Especially, the correlation 
lengths, of which calculation have been one of the major 



difficulty in the string hypothesis, are explicitly evaluated 
in the spin models. For example of the success, we re- 
fer to recent analysis on the quantum-classical crossover 
phenonpip. in the massless XX Z model in "attractive" 
regime. E3l3 

We extend these studies to lattice Fermion systems. 
Our formulation is fully general for the ID Fermion sys- 
tems which are integrable in the sense of the Yang-Baxter 
(YB) equation. As a concrete example, we take the spin- 
less Fermion model with repulsive interactions in gapless 
regime. This simple example already manifests some fun- 
damental differences from the spin models, and yields 
a sound basis for the future studies on more realistic 
Fermion systems such as the Hubbard model. 

As in Refs. |7-20, first one may perform the Jordan- 
Wigner (JW) transformations to the Fermion models and 
further convert the resultant quantum spin models into 
2D classical vertex models. These procedures have been 
successful in studies of the bulk quantities. In evaluating 
correlation lengths, however, this is no longer true. As 
an example, which will be discussed in the main body 
of this paper, let us take Fermion one-particle Green's 
function (CjCi) and its correspondent {crj'a~) in the spin 
model. Obviously they are related, but quite different 
by nonlocal terms due to the JW transformation. At 
zero temperature (T = 0), using the conformal mapping, 
one evaluates the scaling dimensions from the finite size 
corrections to the energy spectra. As the Hamiltonians 
are equivalent through the JW transformations, it is nor- 
mally difficult to discriminate between the energy spec- 
tra of the Fermions and those of the spins. The differ- 
ence lies only in the boundary conditions. Nevertheless, 
even after JW transformation one can explicitly calcu- 
late the correct scaling dimensions only by incorporating 
the proper Fermion statistics at the very last stage (see 
Appendix B). At finite temperature (T > 0), the QTM 
approach gives the correlation function in the spectral 
decomposition form as J2k l^k]"^ i^k / M)^ ■ Here A^ de- 



notes the fc-th largest eigenvalue of the QTM and Ak is 
a certain matrix element. Once the JW transformation 
is performed, it is difficult to trace the difference in the 
statistics in this framework. Then one hardly recognizes 
the difference in the eigenvalues of the QTM between the 
spin models and the Fcrmion models. A simple prescrip- 
tion is not yet found in contrast to the above mentioned 
case at T = 0. 

Let us recall the quantum invepsfi-scattering method 
based on the graded YB relation,tJtJ by which the in- 
tegrability and other algebraic structureS|-af |tiie Fcrmion 
systems have been discussed successfully.E3~E3 

The formulation, however, has a severe problem in ap- 
plying to the finite temperature case. We must treat 
the quantum and the auxiliary spaces on the same foot- 
ing when constructing the QTM. On the contrary, in the 
graded YB relation the quantum space is the Fcrmion 
Fock space, while the auxiliary space is the (graded) vec- 
tor space. 

To overcome these difficulties, we adopt another ap- 
proach to-jtha Fermion systems, which was invented quite 
recentlyrjEa In this method, we consider an _R-operator 
consisting of the Fermion operators alone, together with 
its "super-transposition" . This time both quantum and 
auxiliary spaces are Fermion Fock spaces. Therefore we 
can, for instance, exchange their roles with no difhculty. 
Actually, by careful introduction of the super-trace and 
interchange of it with the normal trace to the partition 
function, we can derive the commuting QTM for the 
Fermion systems. 

The resultant QTM preserves genuine Fermion statis- 
tics. In other words, the selection rule is already built-in 
algebraically. This proper treatment for the statistics re- 
sults in a change of the analytic structure for the QTM. 
In the "physical strip" , the QTM has only one additional 
zero which characterizes "excited free energy" at finite T, 
while in the corresponding spin model there appear two 
such zeros. Consequently, one observes a T-dependent 
oscillating behavior of one-particle Green's function, as 
well as the difference in the correlation length between 
the Fermion model and the corresponding spin model. 
These are smoothly connected to the expected values at 
the CFT fimit, T ^ (see Appendix B). 

This paper is organized as follows. In the next sec- 
tion, we will present the commuting QTM formulation 
of the spinless Fermion model at T > 0. The Fermionic 
i?-operator, together with its "super-transposition" R, 
play fundamental roles. The analytic structure of the 
QTM and the auxiliary functions are discussed in Sec. Ill, 
which leads to the nonlinear integral equations (NLIE) 
characterizing the correlation length. The limit T — > is 
treated analytically at the "half-filling" (rig = 0.5), which 
recovers the prediction from CFT. We also perform nu- 
merical investigations on NLIE and the correlation length 
for one-particle Green's function. To our knowledge, this 
is the first exact computation of the correlation length 
for various interaction strengths, electron filling and for 
wide range of temperatures. In Sec. IV, we comment on 



alternative forms of NLIE derived from different choice 
of the auxiliary functions. They are akin to the standard 
"thermodynamic Bethe ansatz (TBA) equations" from 
the string hypothesis, thus may be of their own interest. 
Details of calculations and supplementary knowledge on 
CFT are summarized in appendices. 



II. COMMUTING QUANTUM TRANSFER 

MATRIX FOR THE SPINLESS FERMION 

MODEL 



In this section we formulate the commuting QTM for 
the spinless Fermion model. The formulation is based on 
the recent developments in the jtudy of the integrability 
of the lattice Fermion systems.L3"L3 The central role is 
played by an operator solution of the YB equation called 
the Fermionic _R-operator. The "transfer matrix" can 
be constructed from the i?-operator, which generates the 
left-shift operator, the Fermionic Hamiltonian and other 
conserved operators. Here and in Sec. II A, we briefiy 
describe the method. 

To extend the method to the finite temperature case 
utilizing the Trotter formula, it is necessary to look for 
another transfer matrix which generates right-shiit oper- 
ator and the Hamiltonian. In Sec. II B, we shall argue 
how to construct the desired transfer matrix by consid- 
ering the super-transposition of the i?-operator. 

Based on these two kinds of the transfer matrices, we 
devise the QTM for the Fermion model in Sec. II C. 
The QTM constitutes a one-parameter commuting fam- 
ily, which is a consequence of the global YB relation. The 
YB relation also enables us to diagonalize the QTM by 
means of the algebraic Bethe ansatz. The free-energy 
and the correlation length are expressed in terms of the 
eigenvalues of the QTM. 



A. Fermionic _R-Operator 

We define the spinless Fermion model by the Hamilto- 
nian 

L 

t r i i 

+2 A (^n, - i) (^n,+i - 1) }, (2.1) 

where c- and cw are the Fermionic creation and annihi- 
lation operators at the j-th site satisfying the canonical 
anti-commutation relations 

{c„cfe} = {c],ct.} = 0, {c],ck}=5jk. (2.2) 



We assume the periodic boundary condition (PBC) on 
the Fermion operators, 



J 



(2.3) 



The parameters t, A are real coupHng constants. In the 
present paper we consider the repulsive critical region 
0<A<1, 0<i and introduce the parametrization 



A :=cos2r/, < 2ry < -. 



(2.4) 



In the subsequent sections, we shall also use the param- 
eter po defined by 



77 

^° ■= 2^- 



(2.5) 



Hereafter we s et t — 1 for simplicity. 

The model (2.1) is exactly solved by the Bethe ansatz 
method. Since the Hamiltonian (E^) preserves the num- 
ber of the particles, we can add the "chemical potential" 
term without breaking the integrability 



't-chcmical • M / ^ 

J = l 



(2.6) 



However we consider only the case ^ = for a while. 

The several physical prope rtie s including the integra- 
bility of the Fermion model ( ^.l| ) has been discussed by 
transforming it into the XX Z model 



H = 



1 ^ 



{aJa^i+aJa^i+A.|a|+,}, (2.7) 



through the JW transformation. However it was rec ently 
discovered that we can treat the Fermion model (2.1) 
only with the Fermion operators. We shall summarize 
the method in what follows. 

First let us consider a two-dimensional Fermion Fock 
space Vj , a basis of which is given by 



|0) 



J' 



|l),:=4lO).' 



c,|0),=0. 



(2.8) 



Define the Fermionic _R-operator acting on the tensor 
product of the Fermion Fock spaces Vj®Vk by 



Ujkiv) := a{v) {-fijnk + (1 - nj){l - Uk)} 
+b{v) {nj{l ~ Uk) + (1 - nj)nk} 



-c{v){c]ck - CiC 



■J'^kh 



(2.9) 



where 



, , sin 77(1; -f 2) , , , sinnv 

a[v) ■■= ^^ , b{v) := ^— -, 

sm Zr] sm Zi] 

A basis of Vj 0Vfc is given by 



c{v) := 1. (2.10) 



|0),®|0)fc:-|0), |l),®|0)fe:=4|0), 

s s ■' 

|0)j®|l)fc:=4,|0), \l)j(^\l)k:=c]cl\0), (2.11) 

s s 

and we can calculate the matrix elements o f (^.91 ) if nec- 
essary. However we keep the operator form ( ^.91) as much 
as possible and avoid the use of the matrix elements, be- 
cause the former is more transparent. ^£^^ -R-operator 
(^.9[) satisfies the following YB equatioi 



7^l2(^t - w)7^l3(u)7^23(w) = 7^23(^')7^l3(^t)7^l2(M - v). 

(2.12) 



The equation (2.12) is an operator identity and one 
should carefully use the anti-commutation relations (2.2) 
to confirm its validity. 

It is one of the fundamental properties of the R- 
operator TZij{v) that TZij(0) = Vij is the permutation op- 
erator for the Fermion operators. 



^j^k^ 



Vjk := (1 - nj){l - rife) - rijUk + c]ck 
Vjk Xj = XkVjk, {xj = Cj or c]). (2-13) 

We can define an analog of the transfer matrix by 

T{v) := Stra {TZaLiv) ■ ■ ■ 7^al(w)} . (2.14) 

Here the super-trace of an arbitrary operator X is defined 
by 



Str,X:=,(0|X|0)a-a(l|^|l)a, 



(2.15) 



where the dual Fermion Fock space is spanned by a{ 
and n,(l| with 



J 



0, 



a(l| 



We also assume 



a(0|0), =a(l|l)a = l. 



(2.16) 



(2.17) 



The super-trace ( [2.15[) corresponds to the PBC for the 
Fermion operators (p. 3D satisfies the property 



StTa{TZaL{v)---TZal{v)} 

= StTa {7^al(^;)7^aL(w) • • • na2iv)} . 



(2.18) 



Hereafter we call (2.14) the transfer matrix for simplicity. 
As in the case with the integrable spin models, the YB 
ensures the commutativity of the transfer 



equation ( 2.12 ) 
matrices (2.14) 



[T{v),T{v')]^0. 



(2.19) 



The expansion of the transfer matrix (2.15) with re- 
spect to the spectral parameter v is given by 



T{v) = T(0) \l + - 



2^ 
sin277 



n+^A]v + Oiv^) 



(2.20) 



which follows from the relationship 
dTZaj{v) 



J^a,j-l 



Av 

-^VajVa.j-l [ Wj-lj + ^A 



(2.21) 



Note that the operator T(0) = Str^ {T'aL ■■■'Pai} is the 
left-shift operator 

T{0)xj = Xj+iT{0), {xj = Cj or c]). (2.22) 



One can easily prove the relation ( 2.22| ) utilizing the 
property of the permutation operator, 



' a,j+ll aj Xj — Xj^il a.j + ll aj, \Xj — Cj Or CJ. 



(2.23) 



B. Super-Transposed Fermionic 7?-Operator 

In this section, we shall consider another transfer ma- 
trix which generates the right-shift operator. For this 
purpose we first define the super-transposition stj for an 
arbitrary operator Xj{v) in the form 

Xj{v) = A{v){l - Uj) + D{v)nj + B{v)cj + C{v)c], 

(2.24) 

by 

Xf'{v) := A{v){l - nj) + D{v)nj + B{v)c] - C{v)cj. 

(2.25) 

Here A{v) and D{v) (B{v) and C{v)) are assumed to be 
Grassmann even (odd) operators. 

Now applying the super -transposition sti to both sides 
of the YB equation ( 2.12 ), we obtain 



nfi (w)7ef2^ {u - ^;)7^23 (v) - 7^23 iv)nf^' {u - v)nfi (u), 

(2.26) 



where we have used a property of the super-transposition 

I). 
(2.27) 



{TZMu)n,i{v)f^ = <^- {v)T^^ (u), {k ^ I). 



Then changing suffixes and spectral parameters as 



1^3, 2^1, 3^2, 

u — > —V, V —> U — V, 



(2.28) 



we get the following new type of the YB equation 

TZi2{u - ^')7^l3(u)7^23(w) = 'R2?.{v)'R.iz{u)'Ri2(u - u), 

(2.29) 

where 



n,u{v):=nt^{-v) 

= a{-v) {-n-jUk + (1 - nj){l - n^)} 
+b{-v) {nj{l - Hk) + (1 - nj)nk} 
^c{-v){c]cl + c,Ck). (2.30) 

Although the new i?-operator TZjk{v) is not symmetric 
{TZjk{v) 7^ TZkj{v)), it is still possible to prove the rela- 
tion 

7^l2(^t - u)7^l3(u)7^23(w) = Tl23{v)Ui3{u)ni2{u - v). 

(2.31) 

Using TZaj (v) 1 we define another transfer matrix by 

fiv) := Str, [UaLiv) ■ ■ -Tlaiiv)] • (2.32) 

Then the commutative properti es of the tr ansfe r matrices 
follow from the YB equations ( ^.29 ) and (2.31), 



0. 



(2.33) 



T{v),T{v')\ = [T{v),T{v') 
The following remarkable relations hold 

Va]Va,j-lXj = Xj-iVajVa.j^l, {Xj = Cj OX c]), (2.34) 

where 



V,k:=njk{0) 



.tJ 



(1 - nj)(l - Uk) - UjUk - {cWk + CjCk). (2.35) 



Using the relations (2.34), one can confirm that the op- 
erator r(0) provides the right-shift operator, i.e.. 



T{Q)xj — Xj^iT{0), {xj = Cj or c'-). 
In other words, T(0) is the inverse of T(0) 

r(o)r(o) = 1. 

Furthermore, from the relationship 



(2.36) 



(2.37) 



Va 



J+1" 



dv 



v=0 



--^7'..,+i7'.,(^H,+i, + 1a 



(2.38) 



the expansion of the transfer matrix T{v) with respect 
to the spectral parameter v is given by 



^(")=^('^)^^-^(^+I^)" + ^("') 



(2.39) 



C. Commuting Quantum Transfer Matrix 



Th e expansions ( 2.2C ) and ( ^.39 ) with the relation 
(^.37) are combined into a formula 



T{u)T{-u) = 1 + ^i|- (h + jA\ u + Oiu"). (2.40) 

This facilitates the investigation the finite tem perature 
properties of the spinless Fermion model (2.1) via the 
Trotter formula, 

exp (-13 (n + ^A j j = ^Ihn^ {T{uN)f{-UN)' ''''^ 



UN 



[3 sin 2ri 
2r]N ■ 



(2.41) 



Here an (even) integer N called the Trotter number, rep- 
resents the number of sites in the fictitious Trotter direc- 
tion and (3 is the inverse temperature (3 = 1/T. 
The free energy per site, for instance, is given by 



/ = - lim lim — - InTr ( T(uAr)T(-MAr) 

L^oo AT— tc» Lf3 V 



N/2 1 

- -A. 
4 

(2.42) 



However, as is the case with the corresponding spin 
model, the eigenvalues of T(un)T(—un) are infinitely 
degenerate in the limit N —>■ go. 

Therefore it is a formidable task to take the trace in 
this limit. To avoid this difficulty, we transform the term 



Tr 



N/2 



(t{un)T{-un)) ' in (|2.42D as follows: 

/ ~ \N/2 

r \T{un)T{-un)] 



N/2 



Tr W Sti 



a2m,0-2m-l 



'^a2^,L{uN) ■ ■ •^a2„.,l('«Af) 



'>^T^a2^-lX{~UN) ■ ■ -Tiai^-iA-UN) , 

Q2™j(^*Jv)^a2™-i,i(-'"w)- (2-43) 



L N/2 

j = l m-- 



We now introduce a fundamental object in the present 
approach called the quantum transfer matrix (QTM) 

TQTMiuN,v) ■.^TrjTj{uN,v), (2.44) 

where the monodromy operator Tj{un, v) is defined by 

N/2 
Tj{uN,v) := Y[ '^a2^.]{v + UN)'R.a2^-i.jiv - Un)- 

(2.45) 



Using the YB equations (2.12) and (2.29), we can show 
that the monodromy operator satisfies the global YB re- 
lation 



TZ2i{v -v')Ti{uN,v)T2{uN,v') 

= T2{uN,v')Ti{uN,v)n2i{v - v'). (2.46) 

Accordingly the QTM constitutes a commuting family 

[TQTMiuN,v),TQTMiuN,v')] =0. (2.47) 

We re mark that the trace in the definition of the QTM 
( ^.44 ) implies the anti-periodic boundary condition for 
the Fermion operators in the Trotter direction,E3 i.e.. 



'^OW+l 



ct.,,,=-4. (2.48) 



'-ai, "-aw + i 



The free energy per site ( ^.42 ) is then represented in 
terms of the QTM as 

/ = - lim lim -LinStr(TQTM(Mw,0)^)-iA. 

(2.49) 

Since the two limits in ( 2.49| ) are exchangeable J3'H we 
take the limit L —f oo first. Because there is a finite gap 
between the first and the second largest eigenvalue of the 
QTM for finite temperature, we can write 

/ = -- lim InAi-^A. (2.50) 

•' /3 N^co 4 ^ ^ 

where Ai is the first largest eigenvalue of the QTM 
TQTMiuN^Q)- From now on A^ denotes the k-th largest 
eigenvalue of the QTM. The correlation length ^ of the 
correlation function (CjCk) can also be represented in 
terms of the first and the second largest eigenvalues A2 
as 



C 



lim In 

N^oo 



A2 

Ai 



(2.51) 



In this way the calculation of certain thermal quantities 
reduces to the evaluation of the eigenvalues of the QTM 
in the Trotter limit {N -^ 00). 

For TV finite, it is possible to diagonalizCr-the QTM 
( p.44| ) by means of the algebraic Bethe ansatzEa (see Ap- 
pendix A). The eigenvalue is then given by 



A{x) ^ Xi{x) + X2ix), 
Xi{x) :— 0+(a;)0_(a; — 2i) 



Q{x + 2i) a^l^ 
Q{x) 



A2(a:) := (-1)^/^+^-0, (x)0+ (x + 2») ^(" ^'K ^^^/^ 



(2.52) 



where 



(/)±(x) := 



f smhri{x ztiuN)\ 



sin 277 



Q{x) := 1 [ sinh77(x — Xj). 



(2.53) 



Here we have changed the spectral parameter from v to 
X defined hy v = ix for later convenience. Note that we 
have also included the contribution from the chemical 



potential term (2^) in the expression ( 2.52 ). 

The associated Bethe ansatz equation (BAE) is given 
by 



(/)-)- (x)0„ (a; — 2i) 
(j}-{x)(l)+{x + 2i) 



N/2 



/3m TT ^{^j 



n 

k=l 



2i) 



Q{xj + 2i) ' 



(2.54) 



Compared with the XXZ mod el, we observe an extra 
factor (-1)^/2+^^0 in (|^ and (g^) which reflects the 
Fermionic nature of the present syste m. I n parti cular , 
if N/2 + TVo = 1 (mod 2), the Eqs. ( |l5^ ) and ( |2^ ) 
are clearly different from the corresponding ones for the 
XXZ model. Actually the second largest eigenvalue lies 
in the sector A'o = N/2 — 1, while the first largest one 
is in the sector N^ = N/2. Therefore the correlation 



length ^ ( 2.51 ) exhib its the manifest difference bet ween 
the Fermion system (p.l[) and the spin system (2.7). 



III. 



NLIE AND THE EXACT ENUMERATION OF 
CORRELATION LENGTH 



A. Analyticities of Auxiliary Functions and NLIE 

In order to proceed further, one needs to clarify the an- 
alytic property of the QTM. For this purpose, we perform 
numerical investigations by fixing the Trotter number N 
finite. 

First we give the description for the largest eigenvalue 
sector, which is naturally identical to the corresponding 
XXZ model. There are N^ = N/2 BAE roots. Only at 
the "half-filling" , they distribute exactly on the real axis 
symmetrically with respect to a; = 0, while for the gen- 
eral particle density they bend in the complex x plane. 
The QTM has N zeros in 3a; e [-2po,2po]: N/2 zeros 
locate on the smooth curve 3a; ^ 2, and the other N/2 
zeros are on the curve 5a; ~ —2. Thus there is a strip 
Qx £ [—1,1] where the QTM is analytic and nonzero. 
We call this "physical strip" . 

Next consider the excited state relevant to the second 
largest eigenvalue. In contrast to the XXZ model, we 
find that two complex eigenvalues are degenerate in mag- 
nitude. Both of them are characterized by iVc = N/2 — 1 
BAE roots located on a smooth curve near the real axis. 
The distribution of the BAE roots for the one and that 
for the other are symmetric with respect to the imagi- 
nary axis. As to the zeros of the QTM, N — 2 zeros are 
on the smooth curves 3a; ~ ±2. 

The locations of the two "missing zeros" are vital in 
the evaluation of the excited states. For the XXZ model, 
both of them enter into the physical strip. Especially, 



with vanishing external field h, they are on the real axis 
and are symmetric with respect to the imaginary axis. 
With the increase of h, they are away from the real axis, 
but still stay in the physical strip preserving the symme- 
try. 

We find a different situation for the Fermion model. 
At the half-filling, corresponding to h = in the XXZ 
model, one of them is located at ^o on the real axis, while 
the other is at 0q + ipo, 9q ^ 6q. Namely only one zero 
appears in the physical strip. Away from the half-filling, 
the zero in the physical strip (we call it 9) moves upward 
while the other {9') moves downward. Nevertheless, we 
find that 9 remains in the physical strip while 9' never 
comes in. From now on we consider the case ^9 > 
{di9' > 0). Then the trajectories of 9', for example, are 
depicted in FIG. 0. 




Re 



FIG. 1. The trajectories of the additional zero 6' are de- 
picted in the case po = 3, A'^ = 100. With the decrease of T, 
8' moves downward, whereas it never comes into the physical 
strip. 

We assume all these features are valid in the Trotter 
limit N —> oo. Then a set of nonlinear integral equa- 
tions (NLIE) can be derived as in the case of the XXZ 
model.Ej We define auxiliary functions 



a(a;) 
a(x) 



Xi{x + i ~ iji) 
\2{x + i ~ iji)'' 
\2{x -i -f 472) 
\i{x -i + 1^2)' 



a(x) := l + a(a;). 



2l(a;) := H-a(x). (3.1) 



where 71 , 72 are small positive quantities introduced for 
the convenience in numerical calculations. Note that 
these functions have asymptotic values 



a{x) 
a{x) 



exp((-7r + Arj)i + (3^) 


for x - 


-> -CX), 


exp((7r -477)1 + 13^) 


for X - 


-^ 00, 



exp((7r — Arjji — (3^i) for x 
exp((— TT + Arj)i — Pfi) for x 



(3.2a) 
(3.2b) 



Immediately seen from the above analyticity argument, 
o(a;),2t(a;) (a(x),2l(a;)) are Analytic, NonZero and have 
Constant asymptotic values (ANZC) in a certain strip 
in the lower (upper) half plane including real axis. The 
above definitions, together with the knowledge of zeros 
for A(a;), fix the NLIE among these auxiliary functions. 
We defer the detail derivation to Appendix C. The re- 
sultant expressions allow for taking the Trotter limit an- 
alytically. Thereby one arrives at the final expressions 
totally independent of fictitious parameter N, 



In a{x) = — 



7r/3 sin 2r/ 



F*ln2l(a;) 



477 cosh -K^ — 171) 

F * ln2t(a; + 2i - ^(71 + 72)) 



In a{x) = — 



■n(3 sin 2rj 



P^Po 



2(po-l)' 
F*lna(x) 

72)) 
2TriT{x ~ e - i{l - ^2)) 



4?7Cosh §ix + 472) 



l3fipQ 



Hpo-1)' 



(3.3) 



where 
A * B{x) 

F{x) := 



1 
2^ 



ZTT 



A{x - y)B{y)dy, 

sinh(po — 2)k 
2 cosh fc sinh(po — l)fc 
sinh(po — 2)k 



-ikx 



dk. 



-ikx 



dk. (3.4) 



, 2fc cosh fc sinh(po ~ 1)^ 

The location of zero 9 satisfies a subsidiary condition, 

a(6l-i + i7i) = -1. (3.5) 

Taking the Trotter limit A'^ — > 00 after setting a; = in 
( pl3D we derive the "first excited free energy " per site 
/2 is 



/2 



P 



lnA2(0) 



-A 



£0 
1 



fi 



K * In2l(i7i) - -K * In2l(-i72 



a 



- — In tanh — 

(3 4 



'2(3 



(3.6) 



where eg is the ground state energy defined in (C2C) and 

f 



K{x) 



4 cosh ^ ' 



(3.7) 



Together with the NLIE for the largest eigenvalue, 
summarized in Appendix C, these relations characterize 
the correlation length ^ of one-particle Green's function 



, Cj) at T > completely (see Eq. ( ^.5l| )) 



We remark that in derivations of above relations one 
does not need precise information like roots distributions 



of the BAE. Only ANZC properties of the QTM and the 
auxiliary functions are sufficient. Thus the structure is 
rather robust, and permits to introduce small free pa- 
rameters 71 and 72. 

In the next two subsections, we present analytical and 
numerical studies on these equations and the correlation 
length of one-particle Green's function, which are main 
results in this paper. 



B. Low temperature property of NLIE (/i = 0) 

We study the low temperature behavior forpUae half- 
filling case /J, = utilizing the Dilogarithm tricked, which 
enables us to obtain the first low temperature correc- 
tion without solving NLIE. As in the case of the largest 
eigenvalue sector, |o(a;)| and |a(a;)| exhibit a crossover 
behavior, 



where 



|o(a;)|,|a(x)| < 1 for |x| < JC, 
\a(x)\, |a(x)| -- 1 for |x| > K, 



^^^2^^.^sin(2,) 



2r/ 



(3.8) 



(3.9) 



Thus one carefully takes into account of contributions 
near "Fermi-surfaces" ±/C. For this purpose, we intro- 
duce following shifted variables and scaling functions. 



la±{x) := Ino ±—x ± /C 1 , 



la±{x) := Ina ( ±—x ± /C ) , 
^:=^(0-/C). 



(3.10) 



and similarly for capital functions 21, 21, A± and A±. In 
T — > 0, they satisfy truncated equations, 

la+{x) = -e-^+7'Ti _^ p^ ^ lA+{x) - F2 * lA+{x) 



+2TiiT\ -{x - 6) + i{l - ^i) 



la+{x) 






(3.11a) 



Fi *lA+{x) - F2*lA+{x) 



-2mT\ -{x-e) -i(l-72) 



(3.11b) 



la^{x) = -e-"^-2'7i + Fi ^ lA-{x) ~F2 'f lA-{x) 

+27riJc-(„oo), (3.11c) 

la-{x) = -e-^+-5»T2 ^p^^ lA_{x) - F2 * lA^ix) 

-27riJ^(-oo), (3. lid) 

where 



F2ix):= -F (-X + 2i- i(7i + 72) 

TT \TT 



(3.12) 



and Fi,F2 are their complex conjugate. In thi s hmit, 
the finite T correction part, lnAf,i(x) (see ( |C14b )) reads 



InAfn a;) - x^+ 2/3 ■ o ' 
2 Tripsin 2?/ 



~27re2 



x-e 



+e" 



(3.13) 



T hank s to the subsidiary condition for the additional zero 
9 (3^), we have 



TT 



2i 



TT 



F{-{z 

TT 



+ i(l-7i) M+(z)dz 



— 00 

F('^(z-0)-i(l-72) )H+(z)dz). (3.14) 

— 00 



TT 



For further simplification, we define D± by, 

/OC 1 7 

(M±(x)-;-/a±(a;) + Z]4±(a:)^-Za±(x) 
_oo V ax da; 

— Za±(x)-— /A±(a;) — la±{x)-—lA±{x) J da; 



'^i^"") /ln(l + a) In a 



a± ( — 00) 



1 + a 



da 



+ 



a± (00) 



a-|- ( — cxd) 



/ln(l + a) Ina 
\ a 1 + a 



da. 



(3.15) 



Obviously, they are equal to special values of Roger's 
Dilogarithm C, 



D+=2C 



-2£ 



a±(oo) 
1 + a±(oo) 

a±(— 00) 
1 + a±(-oo) 



2£ 

-2/: 



£(x): 



1 



dy 



ln(l - y) 



a±ico) 
1 + a±(oo) 
a±(— cxd) 

1 + a±(— 00)/ ' 
Iny 



1-2/ 



(3.16) 



We then apply the dilogarithm trick to ( 3.11a )-( 3.11d ). 
For example, we take the first two equations. After 
differentiating, we multiply them by M+(x),ZA+(x) re- 
spectively and take the summatio n. W e call re sul- 
tant equalit}^ (A). Next multiply ( ^.llaD , ( |3.11bD by 
(M_|_(x))', (/yl+(x))', respectively and take the summa- 
tion. Let us call the outcome as (B). Finally we subtract 
(B) from (A) and integrate over x. The Ihs of the equality 
is nothing but Z5+. Remarkably in the rhs, most compli- 
cated terms like 



M+ (x)^^^^^^-^lA+ {y)dxdy 
dx 

, ^ / N ^ / N dlA+iy) , , 

lA+(x)F2{x - y) P^' dxdy, 
dy 



(3.17) 



and 



dlA+{x) 
dx 



F2{x - y)lA+{y)dxdy, 



(3.18) 



cancel with each other. After rearrangement we obtain, 

D+ + 2tiiT{oo) In i±^ 
A+{oo) 

2e-y {ei''^HA+{y) + e-i''''HA+{y)) dy 



— OQ 
00 



+8i 



— 00 
00 



F[-{x-e)+i{l- 71) lA+{x)dx 



F [ -{x - 9) - i{l - 72) lA+{x)dx, (3.19) 



where a +(~oo) — a^{—oo) = is used. Similarly, from 
( |3.11c| ), ( ^.lld| ), we have 



A-{oo) 
2e-^ {e-^''^HA-{x) + e^'''^M_(x)) dx. 



(3.20) 



Applying ( |3l9| ), ( pO| ), together with ( |3l^ to (|1|) 
2rj 



InAfn(x) r 



2*^ 27r2/3sin277 



An^ + D+ + 2T:iT{oo) In = 



A+(oo) 



-e^T^ ( £>_ +27riJ^(-oo)ln= ^ ' 



A+{oo), 

)}■ (3.21) 



v4_(oo) 
Now that the asymptotic values are easily found 

a+(oo)=a_(oo) = e('^-'"')\ 
a_(oo)=a+(oo) = e(-''+'"')\ 



we can explicitly evaluate (3.21) at x = 



TT TT I 1 a 

6PvF PvF \a 4 



TT 

—i. 
2 



(3.22) 



(3.23) 



where C{x) + C{1 — x) — tt'^ /6 is also applied. Here a 
is introduc ed in (B6) and the Fermi velocity v-p is also 
derived in (B13) for rie — 0.5. The first term is identical 
to the largest eigenvalue sector, and it reproduces con- 
formal anomaly term with c == 1. Comparing them, one 
concludes 



Ai 



(3.24) 



where fcp denotes the "Fermi momentum". Note that 
fcp = 7r/2 in the half-fiUing case. Consequently the in- 
verse correlation length is given as 



r' = 



wp 



(3.25) 



These are nothing but the expected results from CFT 
(see (B14)). This fact represents the consistency of both 
our result and validity of CFT mapping in the finite tem- 
perature problem at low temperatures. 



C. Numerical Analyses on NLIE 

Having verified consistency at the specific limits, we 
now perform numerical analyses on the NLIE for a wide 
range of temperatures, electron fillings and interaction 
strengths. 

To keep the electron filling constant, we adopt the tem- 
perature dependent chemical potential which are deter- 
mined by the curve. 



dT 



^^ll^l ^0 



dT \dfi 



(3.26) 



The NLIE are numerically solved by the iteration 
method. In each iteration steps, convolution parts are 
treated by the Fast Fourier Transformation (EFT). As 
a technical remark, we call an attention to proper re- 
scaling of auxiliary functions for the EFT; one needs to 
modify the integrands such that th ese a sympt otic v alues 
vanish. From the asymptotics in ( [3.2a ) and (3.2b), we 
introduce 



"Six) 



I 2t(x)/2l(oo) for X > 0, 
I2t(x)/a(-cx)) forx<0, 



(3.27) 



and similarly for others. We also rewrite NLIE in terms 
of 23 (x), which has now zero asymptotic values. For ex- 
ample. 



In a(a:) 



nf3sm{2ri) 



4rycosh ^{x — 171) 



F*ln«8(a;) 



-F * In »(a; + 2i- i{-ji + 72)) + T{x) In 

-T{x + 2i- j(7i + 72)) In =-^ 
+2mT{x -e + i{l- 71)) + l3^l. 



a(oo) 

2l(-oo) 



(3.28) 



In addition, one must be careful in the branch cuts of 
the logarithms. In the above, ln2l(oo)/2t(— 00) and so on 
must be understood as 

ln5^^1nr-""'(^^/^-^''^^^l+(vr-4,).. 



"a(-oo) 



sinh(/3/L(/2 -I- 2iri) 



Under these arrangements, the iteration method works 
in a stable manner. 

We plot the temperature dependence of the correla- 
tion length (T in FIG. g for various fillings keeping the 
interaction strength constant A = cos(7r/6). 



0.421 ^-' 
0.398 t - 






(3.29) 



FIG. 2. The temperature dependence of the correlation 
length ^ for po=6- 



The extrapolated values T — > agree with the predic- 
tions from CFT within few percents even far away from 
"half- filling" (rio — 0.5). The curves are going down 
gradually with the decrease of electron density Uq. As 
further information, chemical potential ^{T) determined 
by Eq. ( |3.26| ) and the location of the additional zero 9 are 
depicted in FIG. ^ and FIG. ^, respectively. The zero 9 
moves on a smooth curve and its curvature increases with 
the decrease of ric- In fact, we find that it moves to 9 = i 
when TT-cT -^ 0. (See also the analytic argument for 
non-interacting Fermion case in FIG. 10 for /i = 1.0.) We 
also calcu late the "Fermi momentum" fcp — 91nA2/Ai 
(cf. Eq. (3.24)). (Here the inverse period of oscillatory 
behavior at arbitrary T is referred to as k-p as in the case 
of T — 0.) The figure clearly shows the temperature 
dependency of k-p. In the low temperature limit T — > 
0, it converge to the expected value, fcp = J^eTr, which 
indicates the significance of the Fermi surface for oiife. 
particle excitations in the Luttinger liquid at T = O.cll 
With the increase of T, the auxiliar y fu nctions cease to 
exhibit a sharp crossover behavior ( |3.§| ), which roughly 
corresponds to broadening of the Fermi distribution at 
T > 0. The particle excitations are enhanced within the 
wide range near the Fermi surface, which yield the shift 
of kp . We remark that such T dependent oscillatory be- 
havior has been reported for the longitudinal-correlation 
function of ferromagnetic Heisenberg model.a Although 
the physical origins are different for these two cases, the 
explicit determination of T dependency is important. 




FIG. ^, presents the temperature dependence of the 
correlation length for various interaction strengths for 
fixed Tie. Naturally in the limit, nc,T — > 0, ^T does not 
depend significantly on the interaction strength; it merely 
behaves as ^T ~ vp/ir ~ Uc (see Appendix B). This 
fact is typical for non-interacting cases. Although our 
model inherits strong correlations, FIG. g indicates that 
Tic — 0.1 is already well described by "non- interacting 
approximation" and also shows this approximation is ap- 
plicable in the wide range of T. On the other hand, data 
for Tic = 0.4 show strong dependency on A, therefore 
it belongs to proper "interacting class" (see FIG. 0). It 
seems these crossover occurs near ric ^ 0.25 but it is 
not yet conclusive. We hope to clarify this in a future 
communication. 



FIG. 3. The temperature dependence of the chemical po- 
tential fj, for po~6- 
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FIG. 6. The temperature dependence of the correlation 
length for no =0.1. 



Re 



FIG. 4. The trajectory of the additional zero 9 inside the 
physical strip. 





FIG. 5. The temperature dependence of the Fermi momen- 
tum fcp for po~6- 



FIG. 7. The temperature dependence of the correlation 
length for 11^—0.4. 

Finally we plot the correlation length of trans- 
verse spin-spin correlation {a^aj) without external field 
(FIG. 0) for comparison with rie = 0.5 of spinless Fermion 
models (FIG. y). Besides the difference between their 
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limiting values at T ^ 0, one clearly sees the difference 
in the dependence of ^T on T. 




FIG. 8. The correlation length for {o'^a^ ) of the corre- 
sponding XXZ model with zero magnetic field. 







Po=3 - 
Po=4 - 
Po=5 - 
Pq=6 
analytic calculation 



FIG. 9. The temperature dependence of the correlation 
length at half-filling. 



IV. SUMMARY AND DISCUSSION 

We have proposed the QTM approach to the inte- 
grable lattice Ferniion systems at any finite temper- 
atures. The Fermionic i?-operator, together with its 
super-transposition R, where the Fermion statistics is 
embedded naturally, play the crucial role in this ap- 
proach. Consequently, we have observed the significant 
difference between the Fermion model and that of the 
spin model. In principle, we can apply this approach to 
any integrable ID Fermion systems. The application to 
the Hubbard model is under progress. 

Here we c omm ent on the "attractive regime" t > 0, 
A < in (2.1), which we have not been concerned 



with in this paper. In the XXZ model without exter- 
nal magnetic field, one may recall the remarkable differ- 
ence between the repulsive (anti-ferxopLagnetic) case and 
the attractive (ferromagnetic) one.c3o In the repulsive 
regime, the eigenvalues related to the correlation (cr^cr~) 
or (a^af) is characterized by two real additional zeros 
which are symmetric with respect to the imaginary axis. 
This symmetry is never broken at any temperatures. On 
the other hand in the attractive regime, "level crossing" 
occurs successively. One may attribute it to the change 
of the distribution patterns of the additional zeros. It 
will be interesting to see if similar phenomena happens 
for the spinless Fermion model in the attractive regime. 

Finally we refer to another formulation of NLIE de- 
rived from the different choice of the auxiliary func- 
tions. The NLIE have a close connection with the "TBA" 
or "excited states TBA" equations from the standard 
"string hypothesis" . 

The idea is as follows. First we embed the QTM itself 
into a more general family called T-functions and explore 
functional relations among them (T-system). Then we 
define the ^-functions by certain ratio of the T-functions 
and also derive functional relations for them (y-system) . 
The analytical properties of these functions leads to the 
NLIE which determine the free energy and the correla- 
tion length. As concerns the largest eigenvalue sector, 
the T-functions coincide with those in Ref. 21. There- 



fore the derived NLIE for the free energy are-jidentical 
to the TBA equations of the XXZ model.elEl In con- 
trast, for the second largest eigenvalue sector we find the 
essential difference between the Fermion model and the 
corresponding spin model. For example, we write explic- 
itly the NLIE (excited state TBA equation) for po = 5, 
/I = as 



In 7^1 (a;) 



5/3 sin ^x 
2 cosh ^ a; 



+ K * ln(l + r]2){x) + 7ri, 



ln?72(a;) = K* ln(l - 77i)(l - r]3){x) 

tanh —{x — 6i) tanh —{x — 62) j + m, 
In 773(2;) = K* ln(l + ri2){l - K'^)ix) 

lnK = K *ln{l -ri2){x) + In ftanh j(x - 6*2)) 



2*' 



(4.1) 



where Oi and 62 are determined from 
.5/3sinf6li 



+ K* ln(l + 772)(6'i +i)-m = 0, 

0. (4.2) 



2sinhf6'i 



The meaning of the functions rjj and the quantities 9j 
are similar to those in Ref. ^. Although the above ex- 
pressions are quite different from those in Sec. HI, the 
numerical result shows a good agreement. The detailed 
derivations of above emiations will be described in a sep- 
arate communication e3 
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APPENDIX A: DIAGONALIZATION OF THE 
QUANTUM TRANSFER MATRIX 



Here we shall diagonalize the QTM (2.44) by means of 
the algebraic Bethe ansatz. 

First let us recall that the monodromy operator ( 2.45| ) 
satisfies the global Yang-Baxter relation 



7^2l(w - v')Ti{uN,v)T2{uN,v') 
= '^2{un,v')Ti{un,v)TZ2i{v - v'). 

Writing the monodromy operator as 



(Al) 



Tj{uN,v) — A{v) (1 — rij) + B{v) Cj 

+Civ) c] + Div) n„ J = 1,2, (A2) 



and substituting them into ( Al), we get the commutation 
relations among the operators A{v), • • ■ , D{v), 



A{v)B{v') = ^^^^^—^B{v')A{v) 



b{v' — v) 
c{v' — v) 
b{v' — v) 



Biv)A{v'), 



D{v)B{v') 



h{v — v') 
B{v)B{v') = B{v')B{v). 



*-";ib(„)i,(./). 



(A3a) 

(A3b) 

(A3c) 



To derive these relations, one should pay attention to the 
fact that B{v) and C{v) anti-commute with the Fermion 
operators c, and cJ. 



The commutation relations (A3a)-(A3c) are quite-sim- 
ilar to the corresponding on es for the XXZ modcLEJ In 
fact the relations ( A3a) and (A3c) are identical. The sec- 
ond relation (A3b), however, is different: there appears 



an overall "minus" sign on the rhs. 
Now we define the reference state by 

N/2 

1^):- n |0)-2™®|1)-.^- 



|l)a. 



■=ct 



_JO)a. 



(A4) 



Then using the relations, 

Ra2^-l,3{v - UN)\l)a2^-l 

= -a{-V + UN)n.j |l)a2„-i 

+h{~V + UN){l -rij) |l)a2„-i +Cj|0)a2„_i, 

(A5) 
and 

^a2„j(w + '"Ar)|0)a2,„ = a(v + Ujv)(l - TLj) |0)a2„ 

+ b{v + Un)jIj \Q)a2m ~Cj|1)q2,„, 

(A6) 
we find that 

A{v)\Vl) = {a{v + unM-v + ujv))^/^ \n), (A7) 

D{v)\n) ^ {~b{v + UN)a{-v + un))'^^^ \n). (A8) 



Hence the state \fl) is an eigenstate of the QTM (2.44) 
with the eigenvalue 



smri{v + UN + 2)sinri{—v + un)\ 
sin^ 2r/ / 

sin 77(1; -I- ujv) sin ?7(— w + Mjv + 2) \ 



Ao(«) 



sin 2-1] 

(A9) 

An eigenstate with Nc "particles" can be constructed by 
multiplying the operators B{vj) to the reference state 



Wo 






(AlO) 



Indeed, usingr-thc standard argument of the alge braic 
Bethe ansatz,Ej we can show that the state (AlC) be- 
comes the eigenstate of the QTM if the spectral param- 
eters Vj fulfill the Bethe ansatz equations 

sin ri{— Vj + un) sin ri{vj + UN + 2) 



sinr]{vj + un) sinr]{—Vj + un + 2) 
_ _^_,f/2+N. fj sin,fa--., + 2) 



The corresponding eigenvalue of the QTM ( 2.44 ) 

Tqtm{un,v)\^) = A(w)|*) (A12) 

is given by 

nri{v + UN + 2)sinri{—v + un)\ 



A{v) 



Wo 

n 



sin 2?7 
sin 77(11 — Vj — 2) 



-'- sin?7(u — Vj) 

,f^-.N/2+N^ f sin 77(7; + un) sin-q{-v + un + 2) \ 
^ ' V sin2 2?7 ) 



N/2 



Ne 

n 



sin 77(7; — Vj + 2) 
L sin 77(7; — Vj) 



(A13) 
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APPENDIX B: T < 1 BEHAVIOR AND 
PREDICTION FROM CFT 



Note that p{x) and £{x) are the density function and the 
dressed energy defined by 



We summarize the known results of the correlation 
function a1— XL-.= and its T ^ 1 behavior predicted 
from CFtMeI 

Let us start with the zero temperature case. One- 
particle Green's function|-S|hows an oscillatory behavior 
due to the Fermi surface,cil 



(c^(a;)c(0)) ^ cos{k-px) / X' 



2A 



(Bl) 



The scaling dimension A is evaluated from the energy 
spectra in the finite size system, 



A 



1 



4Z(/Cf 



;{ANf + Z{JCvf{AD) 



(B2) 



Here Z{JC-p) is the dressed charge and /Cp denotes the 
"Fermi surface" satisfying 



R{x) 



R{x - y)Z{y)dy = 1, 
2 sin Arj 



cosh 2x — cos 4?7 



(B3) 



AD and AiV are (half-) integers constrained by a selec- 
tion rule, AD = AN/2 mod 1. For one-particle Green's 
function, they are given by AD = 1/2 and AA^ = 1. 
Thus the critical exponent r]p is defined as 



VF 



2A 



^(/Cf)' 



z{icFy 



(B4) 



The dressed charge Z{JC-p) is explicitly evaluated for two 
special cases ;Ej 



Z{1Cy) = 
where a is 



1 for Tie = (/Cf = 0), 

yapl for rio = 0.5 (/Cf = cx)). 



a 



TT — 2?7 
Then the critical exponent -q-p is given by 



'7F 



1 for Tic ~ 0, 

l/a + q:/4 for rio = 0.5. 



(B5) 



(B6) 



(B7) 



In the scaling limit where CFT is valid, the correlation 
functions at T <IC 1 are recovered by the replacement 



v-p . ttTx 
smh . 

ttT vf 



(B8) 



in the denominator in (Bl). Here vp denotes the Fermi 
velocity 



1 



Vp 



de{x) 



2'Kp(x) dx 



(B9) 



a;— /Cp 



/ N 1 /■ / N / N sin2ry 

^(^^ + 2^ 1_^^ ^(^ y^P^y^'^y = 7r(cosh2x-cos2^)' 

1 f'^'' 

^7^ J-K.F 



sin 2rj 



cosh 2x — cos 277 



/i. (BIO) 



Thus the long distance behavior of one-particle Green's 
function is given by 



{c\x)c{Q)) - cos(fcFx)a;-"''^l"l^/"^ 



(Bll) 



Consequently the correlation length at T -C 1 is identi- 
fied with 



vp 
irripT 



(B12) 



The Fermi velocity ( |B9| ) is analytically calculated for the 
cases rio <C 1 and n^ = 0.5 



.F = r™^ , forne«l, ^g^3^ 

I TT sin 277/477 for 77c = 0.5. 



Therefore we get the explicit correlation length (B12) for 
these two special cases. 



£T=i""' for Tie < 1, ,g^^. 

^ I sin 277/(477(1/0; + a/4)) for 7io = 0.5. ^ ' 



We have also verified the extrapolations from the NLIE 
agree with the prediction (B12). 

Finally we remark on the spin correlation. The main 
contribution to the transverse correlation function simply 
decays algebraically. 



(cr+(x)CT(0)) - l/x 



2A' 



(B15) 



that has no oscillation term. Here A' takes the identical 



form (B2). However we have to use AiV — 1 and AD = 
this time. The difference in selection rules for these in- 
tegers, which originates from the difference in statistics, 
leads to a conclusion 



At^ A' = 



1 



4Z(/Cf)2 



(B16) 



The corresponding correlation length is given by (B12), 
replacing 77F by 773 = 1/2(Z(/Cf))^- One thus obtain 
different correlation lengths simply according to the se- 
lection rules. 
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APPENDIX C: DERIVATION OF NLIE 

For simplicity in notation we define 

c(a;) := a(a; + «7i), <L{x) := 1 + c{x) , 

c(x) := a(a; - 172), €{x) := 1 +'c{x). 

That is, we forget additional shifts for a moment. 
We identify the analytic strips, 



(CI) 



Q{x) 

0_(x) 
(j)+{x) 



3xe(-2po,0) 
3x e [0, 2pq) 
3xe(-2po,0]. 



(C2) 



The following identities are direct consequence of the 
definitions, 



A(a; 



Q{x 



- €(x)- 

"■ ' Q{x - {2po - l)i) 

X(j)^{x + i)0+(a; - i{2po - 3))e~'^''/^ 

A(x - ^)= (-i)A-/2+A-ee(^)Q(^-(2p°-i)') 

Q{x - i) 
x4>+{x ~ i)(j)^{x + i{2po - 3))e'^^/^ (C3) 

Now we consider the second largest eigenvalue case 
Nc — N/2~l. We are in position to utilize the knowledge 
of zeros of A2(a;). 

Consider the integral, 



dz 



lnA2{z)e'^'-dz 



where C encircles the edges of "square" : [zi , Z2] U [z2 , -^3] U 
[23,^4] U [2:4,2:1] in the counterclockwise manner, where 
2i = —00 — i, Z2 — 00 — i, Z3 — 00 + i, Z4 = — cx) + i. 
There is one zero of A2(x) in the region inside C. Thus 
Cauchy's theorem is applied. 



- lA: 

me 



ke 



— \nA2(x~i)e'''^''-'Ux 
ax 



dx 



liiA2{x + i)e''''^'=+'Ux. (C4) 



One substitutes Eq. (C3) into the above equation and 
derives identities among the Fourier components of log- 
arithmic derivatives of Q, £ and £. Explicitly we have, 



dlQ[k] = - 



= fe(Po-l) 



dl<S:[k] 



pfc(Po + l) 



dl€[k] 



4sinh(po — l)fccoshfc 
eH2po~i)dl(f,_[k] + e''dl(f)+[k] 



2 cosh k 



2 sinh(po — l)fc cosh k 
In the above we adopt a notation 

d\TL€.{x) .f,^ 



dl€.[k] := 



dx 



dx. 



(C5) 



(C6) 



etc as the Fourier component of the logarithmic deriva- 
tives. 



On the other hand, from the definition (CI), we have 



dlc[k] = 2dl4)-[k]e''P" sinh(po - l)k 
-2dl(i,+ [k]e-'^^P°-^^^ smhk 
-2dlQ[k]e-^P'>-^^'' sinh(po - 2)fc (C7) 

and similarly for d/c[fc]. 



One substitutes (Cq) into (C7) to obtain a closed equa- 
tion among the Fourier modes of the auxiliary functions. 
Using the explicit form for (/jj-, we get. 



dh[k\ = ~2TTi 



N sinhw^vfe 
2 cosh k 



din 



sinh(po — 2)k 
2coshfcsinh(po — l)fc 



^^- 6^*= sinh(po - 2)k 



+2Trie 



2 cosh k sinh(j)o — l)k 
k+zke sinh(po - 2)fc 



2 cosh k sinh(po — l)k 



(C8) 



By the inverse transformation and integration over x we 
arrive at NLIE. Note that the integrati on co nstan t is d e- 
termined by the asymptotic values in (3.2a) and (3.2b). 
After introducing the shifts 71,2, one obtains the iden- 
tical NLIE in the main text, except for "driving terms" 
as we have not yet taken the Trotter limit N ^ 00. To 
be precise the driving term for lna(a::) is 



N r^ sinW^^,fc(^_,^^)^^^ 



N 
"2 



k cosh k 



(C9) 



Due to the combination of uat — — f3 sin 2ri/2riN and N 
entering above, the Trotter li mit is carried out analyti- 
cally. Then one ends up with (p^). 

The expression for the eigenvalue is derived in a similar 
way. One first notes the "inversion identity" , 



A2(2: + i)A2ix -i) = -'ip{x)(L{xX{x), 



where 
and 

A2(X) 



ip{x) 



(j)+{x — i)(f)^{x + i) 
(j)+{x + i)(l>-{x - i) ' 

A2(a;) 



tanh ^ {x 



-{x + 2i)(j)^{x-2i)'' 



(CIO) 



(Cll) 



(C12) 



is introduced to exclude the zeros of A2 (x) and to com- 
pensate the divergence of A2(a;) at a; ^ ±00. 

Then the Ihs is ANZC in a strip Sx £ [—1,1] and also 
rhs is ANZC in a narrow strip including the real axis. 
One thus can solve ( CIO ) and we get the expression 



In A2(x) = In Ags(a;) + In Afn(a;), 



(C13) 
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where 



APPENDIX D: FREE FERMION MODEL 



TV 



sinh fcu^r sinh(po - ^)^ „-ikx,i. 



, kcoshkainhpok 

+ \n(l)+{x + 2i)(f)_{x-2i), (C14a) 

In Afn(2;) -.^ K * ln2l(x + 171) + K * ln'^{x - 1^2) 



IT TTZ 

In tanh —(x — 0) . 

4^ ' 2 



{C14b) 



Taking the Trotter hmit N ~> 00 after setting a; = and 
using the identity 



hm ln(/)+(2i)(/)_(-2i) 



/3 



(C15) 



we derive the first excited free energy as ( |3.6[ ). 

Next we consider the largest eigenvalue sector Nc = 
N/2. In this case, the spinless Fermion model shares 
same equations with the XX Z model. Then the follow- 
ing NLIE have been already derived in Ref. [l5|. 



lnao(a;) 



ttP sin(277) 



4?7 cosh 1(2: — 171) 



F*ln2lo(a;) 



'F * ln2lo(a; + 2i- i{ji + 72)) + 



2(po-l)' 



lnao(a;) 



tt/3 sin(277) 



477 cosh I (x + ^72) 



'F * ln2lo(a; ~ 2i + i{-/i + 72)) 



F*ln2lo(a;) 

(3fJ,po 
2(po - 1) 



where auxiliary functions ao etc are defined by similar 
way to Eq. ( |3.1| ). Note that their asymptotic values 
a;| — > 00 are explicitly written as 

Oo(x) = exp(/3^), ao(a;) = exp(-/3^). (C17) 

Through above NLIE, Ai(x) is described as 

In Ai(x) = In Ags(a;) + K * ln2lo(a; + 171) 

+K*ln%{x-ij2)- (C18) 

Taking the Trotter limit iV — > 00, we get the free energy 
per site / as, 

/=-ilnAi(0)-iA 

= eo - -^K * ln2lo(i7i) - -3K * ln2lo(-i72), 

P P 

(C19) 
where eg is the well-known ground state energy per site. 



Co 



R{x) 



sin 2?7 



1 



cosh 2r]x — cos 2ri 4 



A. 



(C20) 



Though we do not analyze (C16), (C19) here, they are 
implicitly used in the evaluation of the correlation length. 



Here we consider the free energy and the correlation 
length for the free Fermion model A = (2?7 — tt/2) in 
( ^.l|) . In this case we have (l)±{x ± 4i] = (—1)^ 4>±(x) 
and Q(x±4i) = {-1)^-Q{x) from ( ]2!53| ). Then ^t^) 
simplifies to 



A(x) = g{x) 



Q{x + 2i) 
Q{x) ■ 



(Dl) 



where 



gix) := <t)-{x - 2i)(l)+{x)e^P>' 

+ {-l)^(t>+{x + 2i))4>-{x)e-^f^^'. 



(D2) 



We can easily show that 

K{x + i)K{x -i) = (-!)"•= (?(x -f- i)Q{x - i). (D3) 

This rhs is a known function, which is a distinct feature 
of the free Fermion model. It is convenient to modify the 
function A{x) as 



A(x) = 



A(x) 



(t)+{x + 2i)<f)^{x-2i)' 



(D4) 



(C16) satisfying 



A(a; -I- i)A(a; - i) 
= (-l)f+^" {i^{x) + ^j{x)-'^ + 2cosh{(3fi)) , 



(D5) 



where "tpix) has been already defined in ( pill ). 

First we consider the free energy characterized by the 
largest eigenvalue Ai(a;). It lies in the sector N^ = N/2. 



il)-.j=N/2 



.(1) 



The Bethe ansatz root {x^ }jZi ; ^^j ^ [^1j1] are 
symmetric w ith respect to the imaginary axis. The func- 
tion q{x) in (D2) has N zeros in Sec G [—2, 2]: N/2 zeros 

{xj}jj,i are in the physical strip Qx E [—1,1] and the 

others {x'^jjl^ are out of the strip. As g{x) has a prop- 
erty 



g{x + 2i)\^^{~l)^g{x)\ 



we have 



'-jiM 



•^i \—L 



.(1) 



2i. 



(D6) 



(D7) 



From the BAE (2.54), {xj } are completely equivalent 



to {xj}. Thereby one can shows from (Dl) that Ai(a::) 
does not possess any zeros in the physical strip. 

Since the function Ai(a;) is ANZC Sx e [—1,1], we 
have 

lnAi(a;)== [K *\n {X + X'^ + 2cosh{(3^))] {x). (D8) 
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Using the relations 



Thus we have correlation length ^ for (clci) as 



lim ipix) 

N~>oo 



exp 



cosh^ 



lim A(0) = lim A(0), 
we obtain the free energy per site / as 
/ = -i lim ln(Ai(0)) 

p N-*oo 
1 



(D9) 



= -/ ln(2cosh(/3cosC) + 2cosh(/3^))dC, (DIO) 

in agreement with Ref. |43| . 

Next we consider the correlation length ^ for (dc;). 

The BAE roots {x^ }jj,i~ relevant to the second largest 

eigenvalue are identical with {x^ },-^i except that the 

largest magnitude one xj^i^ = 9 is absent. Then the 

^2{x) possesses the additional zero 6 in the physical strip. 
In the Trotter limit, 9 is given by 



■ sinh 



/3 



n — ij3^ 



(Dll) 



The corresponding zero 0' through the property (D7) is 



9' = - sinh" 

TT 



/? 



TT + i/3/1 



2i. 



(D12) 



The zero 9 {9') never goes over (never comes into) the 
physical strip (see FIG. |0|). 



In 



tanh I — sinh 



/? 



TT — Z/3/i 



^ifsinh-if^^i^Usinh-i 



vr — z/3/i 
~1~ 



(DM) 



In FIG. Owe plot the results (D14) for some fixed par- 
ticle densities. 







FIG. 11. The temperature dependence of the correlation 
length for the free Femiion model. 




FIG. 10. The trajectories of the additional zeros for the 
free Fermion model. The lower zero 9 never goes over the 
physical strip, while upper one 9' never comes into the strip. 
In the case n = 1.0, both zeros moves to 9,9' = i at the low 
temperature limit. 



Consequently, A2 (x) can be expressed as 

|A2(a;)| = Ai(a;)tanh-(a;-6') 



(D13) 
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